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A liquid hydrogen cooled, tape- wound, axially segmented, high-purity aluminum electromagnet 
has been built and tested to 95 kilogauss with a magnet power consumption of 22 kilowatts. Formulae 
for the axial, radial, and circumferential stress components in this type coil are derived under the 
assumptions of isotropy and homogeneity. Both plane strain and plane stress analyses are given. 
The hoop stress is also derived from the "floating shell" concept using thin shell theory. The for- 
mulae from these analyses and thick-wall cylinder theory are used to calculate the hoop stresses in 
the experimental coil, and the results are presented for comparison. The theoretical burst field oi a 
monolithic cylindrical coil disk is derived. 
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1. Introduction 

Many current research programs require the application of high magnetic fields. Plasma 
containment, heating, and acceleration are examples. Cryogenic electromagnets are rapidly 
becoming the most convenient, economical means of producing these high fields. Purcell and 
Payne [l] 1 have discussed the advantages of cryogenic electromagnets and described in detail 
the coil construction treated herein. 

The maximum field attainable with an electromagnet is limited by the mechanical strength 
of the coil. The field strength obtained with a given coolant system and magnet power supply 
may also be limited by the heat transfer efficiency and heat capacity of the coolant and the con- 
ductor resistivity. Therefore, optimum coil design requires simultaneous consideration of stresses, 
heat transfer, and conductor resistivity. The first of these is of sufficient importance to be treated 
separately in this paper. 

In high -field magnets, the mechanical stresses may burst the coil. The coil may also fail by 
localized yielding which can result in increased coil resistance and electrical shorting. There- 
fore, it is necessary to determine the stress field for a specific magnet construction. The coil 
construction of particular interest here is the stack of tape-wound epoxy-bound cylindrical disks. 
The coil disks are alternately stacked with fiber disks which provide electrical insulation and sector- 
like openings for radial coolant passages. The stack of coil disks and separators is shown on figure 
1. This photograph was taken following electrical and intermediate structural failure of the magnet 
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Figure 1. Stack of magnet coil disks and disk 
separators with top four coil disks removed {after 
electrical and intermediate structural failure of 
the magnet at 95 kG). 
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at 95 kilogauss (kG). The foil conductor is assumed to be tightly wound with alternate layers of 
paper insulation and then vacuum impregnated with epoxy resin to form a monolithic cylindrical 
disk. The paper and epoxy combined are not permitted to exceed 5 percent of the disk volume. 
The epoxy-paper tensile and compressive yield strengths exceed those of the conductor. This 
construction technique yields to the simplified analysis of an isotropic, homogeneous cylindrical 
disk. Although the epoxy-bound laminated construction cannot be perfectly isotropic, the assump- 
tion of isotropy is considered representative of the coil construction. These cylindrical coil disks 
are analyzed herein and formulae are developed for the axial, radial, and circumferential stress 
components. The radial and circumferential stress components are derived under the assumptions 
of both plane stress and plane strain. The hoop stress is also developed from the "floating shell" 
concept which assumes that each turn of the coil disk is an individual thin wall cylinder. This 
type analysis should be appropriate for loosely wound, poorly bonded coil disks. The formulation 
for predicting the burst field of a monolithic cylindrical coil disk is also developed. Stress and 
burst field design curves based on the analyses are presented. 

Magnet stresses have been given considerable attention in the literature. For the most part, 
only the circumferential (hoop) stress component has been examined. This is reasonably justi- 
fiable since the circumferential stress is the largest component. Post and Taylor [2] point out that 
the magnetic pressure is proportional to fi?/87r, where B z is the axial field component. They show 
that the main incremental body forces are proportional to the radial derivative of the magnetic 
pressure. Bitter [3] extends the magnetic pressure concept to a very useful form by assuming the 
coil to be a thick wall cylinder with internal pressure B 2 J$ir (Bo is the field at the geometric center of 
the coil). This calculation scheme is quick and easy to use and, as will be shown later, provides 
a rough estimate of the coil hoop stresses. However, it should be noted that this calculation 
scheme is valid only for the hoop stress (o- c ) component. If a shear failure theory is to be assumed 
for the coil, it is necessary to also know the axial (a z ) and radial (07) stress components. If the 
maximum normal stress failure theory is employed, then o> will suffice. It will be shown herein 
that <t z is not negligible and can cause coil failure or excessive deformation under the disk 
separators. 

Detailed stress analyses for various type magnets have been given by Cockcroft [4], Daniels 
[5], Giauque and Lyon [6], and Wells et al., [7]. None of these are deemed applicable to the 
magnet construction of interest for the following reasons: Cockcroft [4] assumed that the mag- 
netic body forces were balanced at the outer periphery by an external pressure. This may be 
true for coils rigidly reinforced at the outer edge but no reinforcement was used for the coils 
described here. Cockcroft also made the simplifying assumption that the solenoid was a ho- 
mogeneous, isotropic body, but neglected the largest stress component (cr c ) in his analysis. In 
addition, the formulae developed are dependent upon the calculation of the mutual inductance of 
a solid coil and a circle in the end plane. Although tables are given for specified variations in 
coil geometrical parameters, the calculation is laborious. Daniels [5] also utilizes the mutual 
inductance of one turn with the entire coil to derive stress formulae for coils fabricated by stacking 
tape-wound cylindrical disks. Since the disks were not considered to be structurally reinforced, 
the problem is similar to the one under consideration in this paper. Daniels [5] considered the 
problem as two dimensional and ignored stresses in the axial (Z) direction. The windings are 
assumed to be homogeneous and anisotropic because the coil was fabricated with nylon insulation 
between adjacent turns, and the composite was not integrally bonded in any way. It is assumed 
that the coil cannot unwind due to friction between adjacent turns. Daniels [5] includes the all- 
important hoop stress and the conditions of stress compatibility in the analysis. It is concluded 
that each turn of the coil is in equilibrium under its own body forces and that axial compression 
does not affect the magnitude of the stresses in the radial and circumferential directions. Daniels 
[5] also concludes that evaluation of the formulae developed is very laborious and suggests experi- 
mentation with models as an alternative to calculating the stresses. Giauque and Lyon [6] treat 
the problem of helical-wound cylindrical coils. Each helical layer of the solenoid was separated 
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by an electrical insulating material which also defined the cooling fluid annuli. Thus, the coil 
could not provide uniformly distributed resistance to radial forces and was prevented from un- 
winding only by frictional drag or deformation of the conductor. Each helical layer of the coil 
is assumed to be self-supporting until the stresses in the conductor exceed the conductor yield 
strength. Giauque and Lyon [6] also use the mutual inductance method in the development of 
formulae for the axial and circumferential stress components. A simplified expression for the 
calculation of the mutual inductance between current element and coil is given for nonhomoge- 
neous coils. This simplified expression is reported to be useful for the quick evaluation of maxi- 
mum stresses in a design feasibility study. Under the foregoing assumptions Giauque and Lyon 
[6] neglect the radial stress component, assume a maximum normal strain failure theory, utilize 
the mathematical expression for calculating mutual inductance, and derive very useful equations 
for the hoop, axial, and failure stresses in unbonded helical-layer cylindrical coils. Wells et al., 
[7] have extensively examined the distribution of magnetic forces in toroidal solenoids. Wells [8] 
has also developed very useful formulae for the computation of hoop and radial stresses in Bitter 
magnets (current density is inversely proportional to the coil radius). The stress formulae devel- 
oped herein are dependent upon the magnetic flux density (B). The recent data of Brown et al., 
[9, 10J provide numerical values of B for any solenoid. Therefore, stress calculations using the 
derived expressions are simplified. 



2. Analysis 
2 . 1 Nomenclature 

B= vector magnetic flux density 
#0= magnetic flux density at the geometric center of the coil 

C= constant of integration 

E= Young's modulus of elasticity 

F= force generated by magnetic field 

/= current per turn of conductor 

,/ = current density in conductor 

L — half length of cylindrical coil 

m= absolute value of the slope of B z /Bo versus y plot 

N = number of turns in coil 
p m — magnetic pressure 

r— radius at some turn in the coil (radial coordinate) 

n = inside radius of coil 

r = outside radius of coil 
t = wall thickness of a thin cylinder 

u= measure of radial displacement 

V= conductor volume 

y— absolute value of BJBq at r= n(y= 1) 

Z = denotes distance along axis from center of coil (axial coordinate). 

a. Greek 

a= ratio of outer to inner coil radius ( = r /r,) 

P= ratio of coil half length to inner coil radius ( = L/n) 

y = ratio of turn radius to coil inner radius ( = r///) 

e = unit strain in conductor 

= circumferential coordinate 

fi= Poisson's ratio 

cr = unit stress in conductor 
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ct u = yield stress of conductor 
4>= magnetic flux geometry parameter 



I a-l Vl + Vl+fi 2 / J 



b. Subscripts 



c = denotes circumferential direction 
r= denotes radial direction 
z = denotes axial direction 



c. Superscripts 
* = denotes plane strain condition 
A = denotes maximum (peak) value 
— = denotes average value. 

The analysis is performed under the following general assumptions: (1) The tape-wound, 
epoxy-bound axially segmented cylindrical coil is considered to be homogeneous and isotropic 
in each disk. (2) The radial and circumferential stresses are derived under the assumptions of 
plane stress and plane strain. The coil is comprised of a stack of relatively thin disks which are 
separated by a disk-shaped insulating material with sector-like openings for coolant flow. There- 
fore, it would be extremely difficult to exactly determine the effect of axial compression on the 
radial and hoop stresses. The assumption of plane stress is judged most applicable since the 
major portion of each disk face is exposed to coolant. Axial stresses due to axial magnetic force 
are transmitted from disk-to-disk via the spokes in the disk separators. The axial compression 
is magnified by the reduction in compression area under the disk separators. (3) The coil is sym- 
metrical with symmetrical flux distribution and uniform current density. (4) Frictional drag 
under the spokes of the disk separators does not resist radial deformation of the conductor. 
(5) Buckling or warping of the disk does not occur. This also means that no bending takes place 
between spokes of the disk separators (due to axial compression). (6) The stress formulae devel- 
oped are applicable only so long as the conductor stress remains below the proportional limit. 
(7) The maximum normal stress failure theory is chosen because the hoop stress component dom- 
inates. The maximum normal strain or maximum shearing stress failure theory is more appli- 
cable to the conductor material immediately below the separator spokes; however, the material 
under the spokes comprises a small portion of the disk volume, and the high compressive stresses 
can be redistributed to the unstressed conductor. 

Consider the cross section of a coil with a magnetic flux linkage as shown in figure 2a. The 
horizontal (radial) component (1) of the vector field induces vertical (axial) forces on the conductor. 
The axial component (2) of the magnetic field induces radial forces on the conductor. The dis- 
tribution of the axial field component across the solenoid radius is shown on figure 2b. The field 
variation along the solenoid axis is shown in figure 2c and the radial field distribution on figure 2d. 
The axial field creates the bursting effect and varies almost linearly across the coil radius. The 
radial field is nearly sinusoidal over the coil radius and increases with distance along the axis (Z) 
for Z < L. Now the incremental force on a conductor segment in a magnetic field is given as 

dV io' (la) 
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Figure 2. (a) Cross section of solenoid showing 
magnetic flux components and associated body 
forces, (b) Distribution of the axial field com- 
ponent over the coil radius, (c) Distribution of 
the axial field component over the coil axis, 
(d) Variation of the radial field component with 
radial and axial position. 



(a) 




0.5 



B r /B 





1 





1 

— r 




Z/L= 


/ / 


,.75 
25 


1 — 


(d) 



The current density for the solenoid in figure 2a is 

NI 



2L(r„ — n) 

Also, the field at the geometric center of the coil may be written, 

D 4tt/V/ . 



(lb) 



(lc) 



Combining eqs (la), (lb), and (lc) we obtain the magnetic body force per unit volume of conductor; 

dF BB 



dV 47T(r — n)<$> 



(2a) 



The radial force component is 



dF r B Z B„ 



dV 4w(ro — n)*' 



(2b) 



and the axial force component is 



dF, B r Bo 



dV 47r(r — r,)$ 



(2c) 
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FIGURE 3. Distribution of magnetic body forces on a coil 
disk at midplane{ZjL = 0). 




Equations (2b) and (2c) are valid for any point in the coil cross section. Equations (2b) and (2c) 
indicate that the radial and axial forces are linearly proportional to the field components. Refer- 
ring to the field distributions in figure 2, we can sketch the loading imposed on a typical disk in 
the coil. The distribution of the magnetic body forces is illustrated on figure 3. 

dF r 
Radial equilibrium of a differential element of conductor with magnetic body force -777 

is satisfied by 



&c -~ cr r — r 



d(T r dF r 



dr 



dV 



(3) 



See Timoshenko and Goodier [11]. 



2.2. Plane Stress 



For the condition of plane stress a z = 0. Hooke's law for the condition of plane stress may 
be written [12] in terms of stress and displacement as follows: 



(T c 



dr 



(1- M 2 )lr ** 
\ du 



du 
~dr 

u 



"(l-/Lt 2 )Ur 

Combining eqs (3) and (4) we obtain the equilibrium equation in terms of displacement 

(Pu ldu u = ( dF r \ (1-/^ 



(4a) 
(4b) 



dr 2 r dr r 2 



dV 



(5) 



Equation (2b) must now be evaluated in order to solve this differential equation. It will be shown 
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later that B Z = B (— — H provides fair results in computing the stresses at the midplane; how- 

Vo-nJ 

ever, to maintain generality we shall let 

B z = B [y+m(\-y)l (6a) 

w here y=the value of B z /B at r=n (fig. 2b), and m = absolute value of the slope of B z /B (fig. 2b). 
Combining eqs (2b) and (6a) we obtain 

dFr_ Bfty+m(l-y)] 

dV 47r(ro-r/)4> ; 



Substitution of eq (6b) into eq (5) yields 

d 2 u , 1 du u T s . /i m ^ , 

-^ + ~J -£ = K y+m(l-y), 7a 

dr z r dr r 

where 

K = ~4nEn(a-m (7b) 

The left side of eq (7a) is recognized as Euler's equation and the complementary solution is written 

c 

Ui = Cir-\ — - (8a) 

r 



The particular solution to eq (7a) is 



„,.^jfcta» 



-*}• 



and u = u\-\- w>. The constants of integration C\ and C 2 are evaluated Irom the boundary conditions 
07 = at the inner and outer radii of the disk. Invoking these boundary conditions and combining 
eqs (4b) and (8) we obtain 

_ KQ,r, KQtf 

C* 1 — ^~J -> C-2" 



24 ' z 24 ' 

where Q t = (y^— ) {^(1 - /*) + 8}<//, + 2) + m(7 + 5/x)}, and 
< ?2-(^ T ) (-Ti;) {8y(/x + 2)(l-a)+m[(7 + 5 i Lt)-8a( M + 2)-fa 2 (9 + 3 A t)]}. 
Substituting into eqs (4) we obtain the following formulae; 
ac = 9Ma- 1) 4> { Ql(1 + M) + g ' il ^ + 3my2(1 + 3 ^ ~ 8y(y + ffl) (1 + 2ja)K (9a) 



9677(a-l)<I> 



;;{<?,(l+/x)-C> 2 ^-^+3myV + 3)-8y(y+m)( M +2)}. (9b) 



The stresses are in dynes/cm 2 when Bo is expressed in gauss. To convert to psi multiply B\ 
X 14.5 X 10- 6 . 
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2.3. Plane Strain 

For the case of plane strain e^ = 0. Pursuing the method used in deriving the plane stress 
formulae we obtain for plane strain, 

B l r^ ™(1-2M) 



967r(a- 1)^-1)0 L "' " z y 



{<?* - Qt K - ^ + 8y(y + m) (1 + /x) - 3my 2 (l + 2 M )J, (10a) 



' 967r(a-l)( i Lt-l)0 r 1 "- y 
where 



0T - Q* . 2 + 8y(r + w) (2 - jn) - 3my 2 (3 -2p)\, (10b) 



Va+l / \or— 1/ cr— 1 

+ 3a 2 (2 M - 3)}, and <?* - ^ {8y(2 - /x) (1 - a) + m[(7 - 2/t) - 8a(2 - /x) 

(a 2 -l)(2/x-l) 

-3a 2 (2 M - 3)]}. 

2.4. Simplified Plane Stress 

If the hoop stress at the midplane only is desired for design calculations some simplifying 
assumptions can be made. 

'ro — A 



Set B Z = B 



ro—ri 



which requires thaty= 1 and m = — ZTV n e( I s W- The resultant equations are 

-=<i^{« i+ *» + ^) + ^-^fe)< i+2 4 <»«» 

where ^ - (y^— ) ffc(l - /*) + 8(/x + 2) + (^~^j (7 + 5/*)}, 

and ^ s (_L_) (_i_) (_£_) {8a(^ + 2)-a^ + 7)-3 M -9}. 



2.5. Bitter Method 

Bitter [3] has suggested the use of formulae for the thick-wall cylinder subjected to internal 
pressure BI/8tt. Utilizing this method we obtain for stresses at the coil midplane, 

_B» l + (a/yf 
ac ~87T a»-l ' (12a) 

_B* l-(alyp 
ar ~87r rf-l ' (12b) 
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As mentioned earlier, eq (12b) is not applicable because the boundary condition at the inner radius 
is not satisfied. 



2.6. Floating Shell Method 

The magnet fabricated at NBS admitted annular cracks of light upon completion, indicating 
imperfect epoxy impregnation. Thus, it seems worthwhile to calculate the hoop stresses as if 
each conductor turn is self-supporting and unrestrained, but closed as in a thin wall cylinder. The 
magnetic pressure may be written 



, dF r (dF r \ , 

dpm= M = Kdv) dr - 



(13a) 



Restricting this calculation to the midplane we set 



b. 



=5o(^y 

Vo - n) 



dr = t, dp m = Ap m , and combine eqs (2b) and (13a) with thin shell theory to obtain 

B\y{a-y) 



<Tc = 



47T(a-l) 2 <J> 



(13b) 



This hoop stress reaches a maximum at the radius r=r /2. Then by the floating shell method 
the maximum hoop stress 



i- \ g! ( a V 



(14) 



FIGURE 4. Equilibrium of a coil disk when the entire 
cross section of the disk is stressed to the yield point. 



dFr 




cr y (r -ri)dz 



o-y(r -r,)dz 





r = r-/r, 
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2.7. Burst Strength of Monolithic Disks 

It is assumed that failure occurs when the coil material becomes plastic throughout. Although 
it has not yet been demonstrated, cr z and o> are relatively low compared to o>, and therefore the 
maximum normal stress failure theory is selected. Referring to figure 4, equilibrium is satisfied 
when 



:-£*, 



2o*j/(ro — ri)dz = d¥ r cos 6. 

(15a) 

Choosing unity thickness {dz = 1) and combining eq (6b) and (15a) and performing the integration 
we obtain the theoretical bursting field 

1.73(a-l)WP 
°~ 5~ . (15b) 

(y+ m )(a 2 -l)-y(tf-l) 

where B is in kG and cry in pounds per square inch. 

If for the sake of simplicity we substitute y= 1; m= l/(a— 1) into eq (15b), the bursting field 
at the midplane is 

* i== (« + 2) ■ (15C) 

2.8. Axial Compression 

The axial magnetic body forces are transmitted from disk to disk by the spokes of the disk 
separators. Thus, the forces are additive from the ends to the center so that the disk in the mid- 
plane must support the axial body forces of the entire coil. Axial equilibrium of a differential 
element requires that 

dor z dF z 

1^ = W (16) 

In order to obtain cr z we must evaluate B r (r, z) in eq (2c). It will be shown in the sample calcula- 
tions which follow that the radial field component is approximated by an expression of the form 

^ = Ko + (k-K ){ S mir(^\)}, ( l 7a) 

where K = \ — —7: — — \ for each Z/L of a fixed coil geometry, and K = I 7^-) (peak 

I 2 J \Bq/ g+i 

2 

value of B,/Bo for each Z/L of a fixed coil geometry). K and K — Kq can be linearly related to the 
axial dimension Z by; 

K = m x ZIL, K-Ko=m 2 Z/L, (I7b) 

where m\ is the slope of the K versus Z/L plot, and m 2 is the slope of the (K — Ko) versus Z/L plot 
for a fixed coil geometry. Combining eqs (I7a) and (I7b) we obtain for the radial field, 



t=zh+ m * 8in *(H> (l7c) 
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Combining eqs (2c), (16), (17c) and performing the integration the axial stress is derived, 

m 2 



87rL/-/(a-l)* 



ni\ 



+ m 2 sin 7T H^3j) [+/W + C - (18a) 



Since external forces are not being considered C = and/(r) is obtained from the boundary con- 
dition (cr z )z=L = 0. Performing the indicated algebra we arrive at a solution for the axial com- 
pressive stress 

<7z = "87rL>-l)d>h + m2Sin7r l^TJr (18b) 

This expression is applicable throughout the axial dimension of the coil. The axial stresses are 
largest at Z = 0, and 

{a * ) ™ = -ma-iYt>{ m, + ms sin "(^l)}' (18c) 

a+ 1 
The radius at which (ct z )z=o becomes maximum is y = — - — Assuming that each disk acts as a 

rigid plate so that the axial stress is equally distributed we find the average axial stress (a z ) from 

_ 1 f a 

Vz = ——[ Vzdy . (19a) 



Performing the indicated integration on eqs (18b) and (18c) we obtain 

and 

((r ^ -"87r(«-l)*r i+ ^rj- (190 

In order to obtain the maximum average compressive stress under the spokes of the disk separa- 
tors at the midplane we need only to multiply the results of eq (19c) by the area magnification ratio 
(^diskM spokes)* 

2.9. Sample Calculations 

In order to compare the various formulae and illustrate their use, sample calculations are 
performed for the coil geometry of the NBS magnet. 

a = -=3.67; £ = -=1.87; r, = 1.5in.; r = 5.5in.; L = 2.8in.; /x, = 0.33 (aluminum); ^ = 4500 

psi; = 0.64. To evaluate eqs (9), (10), and (18c) we must determine the values of y, m, mi, and 
m 2 . All calculations will be performed for the midplane (Z = 0) where the stresses are largest. 
All data for the radial and axial field components are taken from the tabulations of Brown et al., 
[9]. Brown and Flax [10] recently presented a graphical method for determining the magnetic 
field components of any finite coil by superposition of four semi-infinite solenoids. The axial field 
component {B z jBo) is plotted for various coil parameters (a, /3) in figure 5. The dotted line indi- 
cates the extrapolated curve chosen to fit the coil parameters of the NBS magnet. Since /3 is 
shown to have little effect on the slope, the selection of the dotted line is based on linear extra- 
polation between the curves for a: = 3 and a = 4. From figure 5 we obtain y— 1.04 and \m\ =0.45. 
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fi = 2, a = 3, Z/L = I 
/&=2,Ct = 4,Z/L =1 
/?= I , a = 4, Z/L --/ 




Figure 6. Variation of end-plane radial field 
with coil geometry. 



7 ' r/r| 

FIGURE 5. Variation of midplane axial field with 
coil geometry. 



On figure 6 is plotted the radial field component for various coil parameters. It is difficult to extend 
these curves to intermediate values of a and /3. The effect of smaller f3 is to increase the magni- 
tude of B r /Bo with little phase shift. Thus we plot on figure 7 the B r IB data for a = 4,/3 = 2, and 
various Z/L. We note that the curves are essentially sinusoidal from y=ltoy = 4 = a. Recalling 
that 






Ko = 



B r 

,6o/y=i 



and 



K 



0/ y=a±l 



we plot Ko and (K — Ko) against Z/L on figure 8. From the Ko plot we obtain the slope mi =0.213, 
and from the (K — K ) curve we obtain the slope m 2 = 0.l5. Substituting these values of mi and 
m,2 into eq (17a) we check on the validity of the assumptions. The B r /Bo is calculated from eq (17a) 
and plotted on figure 7 for comparison. The calculated values plot high for lower values of Z/L; 
however, recall that the effect of a lower /3 (fig. 6) is to increase the B r /Bo. Also, the fit is good 
for higher Z/L where the largest body forces are generated. The various formulae are evaluated 
and plotted for comparison on figure 9. From eqs (15b) and (15c) we calculate the bursting field 
as 93 kG and 84 kG, respectively. Using eq (19c) and an area magnification ratio of 4.4 we obtain 
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□ - FROM Eq (17a ), 

m, = 0.2/3, m 2 = 0. 15 
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FIGURE 8. K and (K — K ) as a function of axial 
position for a fixed coil geometry (a — 4, /3 = 2). 



FIGURE 7. Radial field as a function of axial position 
for a fixed coil geometry (a — 4, j3 = 2). 



the average axial stress (psi) under the separator spokes as —0.843 B\ (Bo in kG). The maximum 
hoop stress [(cr c )y=i from eq(lla)] is plotted for various coil parameters (a, /3) on figure 10. Figure 
11 shows the bursting field [B from eq (15c)] as a function of a and /3. 



2.10. NBS Magnet Failure 

The NBS tape-wound aluminum magnet failed electrically with intermediate structural dam- 
age at 95 kG. The magnet failure is attributed to excessive deformation at the outer radii which 
normally precedes bursting. Therefore, the failure field should not be much different from the 
actual burst field. The failure field agrees well with the 93 kG burst field calculated from eq (15b). 
The actual burst field is expected to exceed the theoretical value due to the frictional restraint of 
the axial compression (which was neglected in the analysis). Eight of the fourteen disks suffered 
circumferential cracks at radii of about 3.75 in. indicating excessive deformation at the outer radii 
and structural inhomogenieties. Some of the cracks were 0.125 in. wide. There is evidence of 
some radial deformation over the entire coil radii. The deformation of individual coil disks varied 
with coil axial position. The disks nearest the coil midplane were exposed to the largest magnetic 
pressure and suffered greater damage. The pronounced failure at an intermediate radius has been 
attributed to an intrinsic weakness in the coil disk construction. Several of the disks were rather 
loosely wound, and consequently some sizeable voids in the outer turns of the disk were epoxy- 
filled. All of the disks that failed possessed these voids from which the circumferential cracks 
propagated. Other disks in the coil had been tightly wound and suffered little or uniform stress 
damage. It would be quite difficult to determine what type structure the loosely wound epoxy- 
bound disk construction represents; however, it is estimated that the stresses should be interme- 
diate to floating shell and plane stress analyses. The hoop stresses for the floating shell (eq (13b)) 
and simplified plane stress (eq (11a)) conditions are superimposed on the conductor yield strength 
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Figure 9. Circumferential, radial and axial stress distributions at the midplane 
for the NBS magnet (a = 3.67, ^8 = 1.87) as computed from the various theories 
discussed. 

on figure 12. The variation in radii at which failure occurred is also indicated on figure 12. The 
conductor yield strength (0.2% offset) was experimentally determined to be 3300 psi at room 
temperature. The yield strength was then inferred to be 4500 psi at 20 °K (the operating temper- 
ature of the magnet). The actual strength of the material throughout the coil was unknown because 
the magnet had been tested at fields up to 70 kG several times prior to the 95 kG test. Thus, the 
material in the inner turns had been strain-hardened. Increased magnet resistance substantiated 
strain-hardening. The cyclic effect of strain on resistance and strength is under study; however, 
it is well known that the yield strength of pure aluminums can be increased three-fold at 20 °K by 
strain-hardening. A three-fold increase in yield strength of the inner turns of the coil disks would 
explain the circumferential cracks at intermediate radii when the magnetic field was increased 
from 70 to 95 kG. The epoxy voids provided the inherent weakness to permit the outer turns to 
expand radially outward from the strain-hardened inner turns. From figure 12 one could only 
conclude that the loosely wound turns of the epoxy-bound disks conformed more closely to the 
floating shell analogy. The shaded area on figure 12 indicates the region where the actual stresses 
for the loosely wound coils may plot. Some of the disks near the midplane were also badly dam- 
aged by axial compression. Extrusion of the conductor into the flow passages varied from 0.001 
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FIGURE 10. Maximum hoop stress as a function of 
the coil parameters, a, /3. 



FIGURE 11. ^arsf ^e/J as a function of the coil 
parameters, a, j8. 



to 0.014 inches depending upon the disk axial position and the apparent aluminum-epoxy com- 
posite strength. Although the disks are separated by 0.025 inch fiber spacers a 0.004 in. feeler 
gage could not be inserted in some of the flow passages upon post-run inspection. 



3. Conclusions 

For quick design computations the method of Bitter, eq (12a), though nonconservative, is by 
far the most convenient. If a more precise analysis of both circumferential and axial stresses is 
desired eqs (9a) and (18c) are recommended. Equation (11a) is more convenient than eq (9a) and 
provides conservative results for hoop stress calculations. 

Equation (13b) is suggested for hoop stress calculations on coils which are loosely wound and 
not integrally bound. The radial stress component is so small that it can be neglected in the coil 
construction described herein. Equation (lib) is recommended for those coil constructions where 
the radial stress is pertinent. Equation (15b) predicts the burst field of the tape-wound coil disks 
with reasonable precision. Equation (15b) and the simplified eq (15c) predict burst fields lower 
than those to be expected in practice. The actual burst field of a coil may be increased by strain- 
hardening the conductor through cyclic testing to the requisite fields. 
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Figure 12. //oop stresses (as computed from the 
most dissimilar theories) and yield strength are 
superimposed on the failure radii of the NBS 
magnet. 
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